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O ■ In this and the companion paper a novel holonomy formulation of so called 

^ ■ Spin Foam models of lattice gauge gravity are explored. After giving a natu- 

r^. ral basis for the space of simplicity constraints we define a universal boundary 

Q-i Hilbert space, on which the imposition of different forms of the simplicity 

^ \ constraints can be studied. We detail under which conditions this Hilbert 

, ^ space can be mapped to a Hilbert space of projected spin networks or an 

ordinary spin network space. 

^ . These considerations allow to derive the general form of the transfer oper- 

Q\ '. ators which generates discrete time evolution. We will describe the transfer 

^ I operators for some current models on the different boundary Hilbert spaces 

t::J- I and highlight the role of the simplicity constraints determining the concrete 

Q\ • form of the time evolution operators. 

o: 

(N. 

7^- 1 Introduction 

>■ 

^ I Together with the companion paper [15] this paper introduces and studies a class of 

■ lattice gauge theories that ocurr as spin foam models of quantum gravity [891 EHl [IDl ESI 

[III El [791 [THl lai Ea ESI EZl [ig [sa lag 

As opposed to previous formulations that stressed the relationship to the loop quantum 
gravity Hilbert space [031 [SlISl [313 IHH [H] , and thus was given in terms of spin networks 
and their geometric interpretation, we will here focus on a presentation that is as close as 
possible to lattice gauge theory. There is some overlap with the formulations explored in 
the context of auxilliary field theories on the group manifold, [ ^ ISO j 19011911 Ii21[88ll92[ 
[H21 [SD IS31 [21] 5 though our perspective is markedly different. Further there are similarities 
to the formulations of spin foam models in connection variables [HSl [HSl [HSl [SZl ISB [ZH ISO] 
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and as more ordinary discretized theories [561 EHl EZl [351 [36], and non-commutative first 
order formulations |20l [21]. In a fortlicoming work [65] we will show how to extract 
geometric meaning directly from the formulation given here. 

The formulation we explore here takes as its starting point the heuristic formulation of 
lattice BF theory on a 2-complex. We then insert simplicity constraints in the holonomy 
picture of this formulation, assuming the most general form common to most simplicity 
constraints. This allows us to define a very natural space of simplicity constraints 
covering almost all models in the literature, the exception being the model due to KKL 
[nSl[HZ|- Combined with the results on the structure of the simplicity constraints of spin 
foam models in [65] this enables us to give natural extensions of the EPRL/FK and BC 
model to arbitrary, including finite, groups that we explored in |15] . 

The key aim of this paper is to explore the consequences of this generic formalism. 
We will give an explicit basis for the space of simplicity constraints in section 12. 3^ and 
discuss which models lead to real partition functions. In section [HI we then explore 
different ways of introducing boundaries and composing 2-complexes. These will lead 
to different notions of boundary Hilbert spaces. Next to the known spin network [9l] 
and projected spin network [3l [72] spaces we introduce a new universal boundary space 
common to all spin foam models built on the same group. We will also show how the 
assumptions we make on the structure of the simplicity constraints translates to the 
operator spin foam formalism, show how the usual models fall into this category and 
give the basis coefficients of the established models in the space of simplicity constraints 
in section [5] 

In section [6] we use the results obtained so far, to derive the general form of the transfer 
operator for any holonomy spin foam model in the different boundary Hilbert spaces. 
Transfer operators arise from a space-time decomposition of the partition function and 
indeed generate the (discrete) time evolution on the given Hilbert space. Thus transfer 
operators can be seen as the discrete time equivalent of Hamiltonian operators. Indeed 
for standard lattice systems Hamiltonians can be derived from transfer operators by a 
limiting procedure. We will shortly touch on the issue of how to take the limit in section 
[HI The general form of the transfer operator allows to illuminate the dynamics of spin 
foams, to highlight the role of the simplicity constraints and to clarify the connections 
between the different boundary Hilbert spaces. 

We end with a discussion and outlook in section [3 

2 Holonomy Formulation of Spin Foam Models 

We will begin by recalling the formulation of spin foam models based on two arbitrary 
finite or Lie groups G and H d G m terms of holonomies on an arbitrary 2-complex. 

2.1 The Data 

We start with an arbitrary, finite, combinatorial 2-complex C. The 2-complex C consists 
of faces, edges and vertices and we write f ^ Cf, e ^ Ce and f G C„ respectively. 
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Each face of the 2-complex comes with a fiducial orientation given by the order of edges 
and vertices around it, as well as a fiducial base vertex. This orientation is unrelated 
to any orientation on a manifold from which the 2-complex might be constructed. The 
combinatorial information, together with the fiducial orientation and base vertex can 
be conveniently encoded by thinking of a face as an ordered set of the vertices and 
edges bordering it, / = (w, e, w', e', . • • v). The notation (a, 6) C / will always mean 
that the ordered set (a, b) exists as an uninterrupted subset in the ordered set /, that 
is / = {v, . . . ,a,b, . . .v). We will similarly write w G e and e G / to denote adjecancy 
relationships. 

As data on these 2-complexes we introduce one G element gev = Qve P^r half edge, 
and one G element gej per neighbouring edge and face, e G / □. In order to have a 
natural composition of group elements this should be read from right to left, that is, 
the group elments Qah and g^yc associated to the ordered sets (a, b) and (6, c) respectively 
compose naturally to gab9hc if we read Qab as going from 6 to a and (7bc as going from c 
to b. 

From these data we define two types of face holonomies as, 

9f = n 9abgbfgbc 

(a,6,c)C/ 
beCe 

9f = Yl 

(a,fe)C/ 

Or, more eliptically, gf = gvege/gev'gv'e'ge' f ■ ■ ■ and gj = gvegev'gv'e' ■ ■ ■ 
2.2 The Partition Function 

We will first consider the partition function for a 2-complexes without boundary, that 
is, we treat all edges and faces as internal. We will see later that most spin foam models 
considered in the literature fall under the following definition: 

Definition 2.1 (Holonomy Spin Foam Model (no boundary)). Let C be a 2-complex 
with orientations and base points on the faces, G a unimodular Lie or discrete group 
with Haar meassure dg, and H G G a subgroup of G. Then given two complex valued 
distributions E and a; on G satisfying 

• E{hg{hy^) = E{g) for all heH, 

• E{g) = E{g~^) and u{g) = uj{g-^), 

• ^{g) = ^{ggg'^) for all g eG, 

^Note that these gef have nothing to do with the holonomy from the mid point of the edge to the 
midpoint of the face that is introduced in the wedge formahsm. 
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we define the Holonomy Spin Foam partition function defined by E and u to be: 

^(^) = / (n^^^^) (n^^^'^) (n^(^^^)) (n^(^/)) • 

If we furtlier liave u}{g) = oj{g) and E{g) = E{g) we say tliat we liave a Real Holonomy 
Spin Foam Model. 

An immediate consequence of the properties of E and u is that Z[C) is independent 
of the orientations and base vertices chosen. 

Note that the partition function as given will usually diverge, even for compact groups. 
For noncompact groups there is also an ambiguity in the normalisation of the Haar 
measure and additional divergences due to gauge orbits. However, if w is a regular 
function and the group is compact the model will be well defined. For finite groups this 
is always the case. 

In the rest of this section we will study some general consequences of this defini- 
tion, arising mostly from the structure of the integrand. Thus we will ignore issues of 
convergence from here on. 

2.3 A basis for the space of theories 

Note that as the conditions on E and u are linear, the space of partition functions given 
two groups H <Z G carries a linear structure. Thus we can parametrise it by giving an 
explicit basis for the space of distributions E and u. 

The distribution is a class function, thus the characters form a basis, and for compact 
groups it can be expanded as 

u{g) = Y,^im{p)u''tip{Dp{g)) , (3) 
p 

with Dp{g) the representation matrix of g in the unitary irrep p. 

We will often set u = 6g, or u'' = 1 in which case the space of theories is simply 
parametrized by the functions E satisfying the conditions of definition 12.11 

E encodes the analogue of the simplicity constraints for the spin foam model at hand. 
The presence of the delta function on the face ensures that the product of group elements 
around the face, gf, is fiat. This is however not the usual holonomy around the face gj, 
but cjf interwoven with the gef- If we force gef = 1 we have gf = gf, and obtain a theory 
of fiat connections. The presence of gef and functions E that allow them to differ from 
the identity thus relaxes the constraints on flatness. 

Thus they exactly play the role of the simplicity constraints in ordinary spin foam 
models. We will see the precise relationship between the simplicity constraints in the 
usual operator spin foam models and the functions E in the next section. We call this 
space of simplicity functions S{G, H). 

We can given an explicit basis for this space by expanding the functions in terms of the 
irreducible unitary representations of G and H, which we denote p and k respectively. 



4 



We have chosen for every irreducible representation some specific reahzation. In the 
space HorriHip, k) we can introduce scalar product by 

Homnip, k) x Homnip, k) 3 (J, /') (/, /')Hom^(p,fc)lfc = I^I' ■ (4) 

For every p we have an antiunitary, group covariant map 

Jp-.p^p (5) 

If p and p are distinct we assume that Jp = J"^ . If p = p then we may assume in 
addition that = Splp (where Sp = ±1). 

Similarly, for every k we have antiunitary map 

Jk-.k^k (6) 

If k and k are distinct we assume that Jj, = J^^. If k = k then we may assume in 
addition that = s^lfc (where Sk = ±1). 

We can thus define the antiunitary (in the scalar product (jl])) map 

Mp^k : Honinip, k), -> Homnip, k), Mp,k{I) = JpUl (7) 
Here J"!" is defined by 



(■,jt.) = (J,.). (8) 
In the case when p = p and k = k 

^p,k = '^kSp lHomH{p,k) (9) 
±1 

Let us define 

D{g)p,k ■■ Homnip, k) Homnip, k) (10) 

by the matrix elements 

{I,b{g)p^j')=tTpDp{g)r . (11) 
We can then expand E as such: 

E{g) = dim(p) tTHomH{p,k) elbp^k{g~^) (12) 

p,k 

where : Homnip, k) — )■ Homnip, k)- ^{G, H) can then be parametrized through e^. 
We will need the following set of useful relations satisfied by D. We see that 



(/, Dp,k{g~')l') = tr PDl{g)l' = W)Wp{g)I = ^^^^ 

= {l',DpMl) = {bp^kig)l,l') 

thus bp^k{g'^) = bp^kigV- 
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Similarly 



{I,Ml,D,-,{g)M,4n) = {M,4l),D^;,ig)M,4l')) 



(14) 



= tiJ, Jt D,{g) 4 J, I' 4 = {I,D,,kI') 
where we used identity valid for any antiunitary J and linear A 

ti J AJ^ = trA (15) 

Thus Dp4g) = Ml,D^-,ig)M,,,. 

The are not completely free, but are restricted by the condition that E(g) = E(g~^). 
This implies that they have to satisfy a set of equations relating the coefficients for 
complex conjugate representations. By definition we have that 



tTe',D,^,{g~') = trei^DpM = tr M;,D,,s(<7)M,,, = tr M,,,e^ Ml,D-^-,{g) (16) 

The condition for the E function read^ 

4 = Mp^ue^Ml, . (18) 

Let us notice that if p ^ p oi k ^ k we can choose bases in such a way that Mp^k acts 
by complex conjugation, then 

el = ef (19) 

The same is possible if SkSp = 1 but not in the case when SkSp = — lH. From now on we 
will assume that SpSk = 1 whenever p = p and k = k. In many models we will consider, 
p and k are indeed isomorphic to p and k, in which case the condition simply says that 

has to be symmetric. 

Furthermore as 



tre'^Dp^.ig-') = tre^pM^ = tie^DpM (20) 
we have a real holonomy spin foam model if 

< = . (21) 
In the cases where we have (11 9p this reads 

< = 4- (22) 



^Matrix elements of Dpj- satisfy 

dp [ dg{h,bp^k{g-^)h){h,bp^k{g)h) = 5p.p.5k,k'dk{h,h){hj2) (17) 



thus the basis form a set of independent functions of g. 
■^If the degeneracy is 1 then always s^Sp = 1 since the restriction of Jp to k is equal to Jfe 



6 



We can write these conditions exphcitly by choosing an ortho normal basis in Homnip, k) 

I{p,k)d (23) 
satisfying conditions from above. In particular this means 

I{p,k)i,I{p,k)d = 5dd'lk (24) 
Using formula E{g) = ^^dim(p) tr e^Dp k{g) and f ll9p the E function can be written 

as 

E{g) = J2 dim(p)e^,,,, tr, {D,{g) I{p, kU\p, k),,) , (25) 

S{G,H) can then be parametrized through e^^^,. The index d can be seen as a 
degeneracy index for the case where the H reducible representation p contains more 
than one copy of the irreducible k. 

3 Boundaries and Hilbert spaces 

We can now introduce a notion of boundary, which will lead us to a new, and more 
general notion of boundary Hilbert space. 

3.1 Boundaries 

We can now introduce and study boundaries into the formalism. To do so we chose an 
arbitrary graph T in C, with edges Eg C Ce and vertices E^, C C„ subsets of the edges 
and vertices of C respectively, as the boundary graph of C. We then write Ee„ for the 
set of pairs v e e in E„ x Eg. 

The partition function is then an element of the space ((?!'"'=" ') by dropping the 
integration over the group elements associated to these pairs: 



Z^iC)[ge,] = 




This definition has the advantage that the inner product of the partition functions cor- 
responds to the gluing along the graph. That is, for two complexes C, C with isomorphic 
boundary graphs we have: 

(WiC),Z^{C')y = Z{CUrC') (27) 

where C Ur C indicates the 2-complex with E in both complexes identified and now 
internal. This follows immediately from definitions. Eor a real spin foam model we of 
course can drop the complex conjugation. Note that the edges and vertices of E become 
edges and vertices of C Ur C. In particular these edges can be bivalent. 
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The integrand in (j2]), given by 



u^idef)] n^(^/))' (28) 

\eCf J \ f 



has the following symmetries: 



9ev ^ QevQv 

Qef he'^Qefhe (29) 

for he & H G G and G G. 

Due to the symmetries of the integrand the partition function actually can be consid- 
ered to live in a smaller subspace of L'^(G^^^^^), that is, 

Zen"^ = L^ (gI^™! /(GI^^I X //l^^l)) , (30) 
with the action of (gv.he) G (G''""' x ifl'""^') on g^^ by left and right multiplication: 

We call this space the universal boundary space for the class of models £{G^ H) and 
write 'Hij^g. 

This should be contrasted with the usual projected spin networks space 

Ze'H'^Psr, = L'{G\^^^/{H\^^^)). (31) 

3.2 The spin network basis 

We will now briefly give the spin network basis for the space 'Hubs- 
It is convenient to start with the basis for the larger space 

L2(GI^^"I)- By the Peter- 
Weyl theorem this is given by the matrix elements of representations, that is 

\p,ij) = ^/d^mJJJ}Dp{g)ij , (32) 
with the dimension factor providing the correct normalisation, 

{p',i',j'\p,hj) = V dim(p) dim(p') j dgDp,{g)^,.,Dp{g)ij 

= ^J dim(p) dim(p') j dgDp,{g^^).,^,Dp{g)ij 
= 5ii/Sjj/5pp> . (33) 
A basis is thus simply given by the tensor product of basis elements 

\pev,ievjev) = H V ^Mpev) D p^^,{9evh^j^^ G ^^(GI^-I) . (34) 
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T-LiTBS is subspace of states in L^(G''"'="') that are invariant under the action of the 
symmetries. In order to give a basis of this subspace it will actually be more convenient 
to use the orientation on the edges to introduce the oriented basis. For this we choose 
an arbitrary orientation for each edge. We will encode this by writing (w, e, v') G F for 
the oriented edge e running from vertex v' to v. 

jev) ■ (35) 

e 

{v',e,v)Gr 

We now want to go to the G invariant subspace of the pev at the vertices, and the H 
invariant subspace at the edges. We begin by implementing the invariance under H. To 
do so we contract the inner indices at each oriented edge with an H covariant operator. 
These are parametrized similarly to the E functions, by the matrices Sk^dd' which can 
be contracted with a basis of H invariant maps between p and p' giving 

En: = Ek,dd'IiP^k)dAmHp'^k)l,^^-,. 

k,d,d' ,m 

Contracting these on the indices i^v in the middle of the edge gives 



\pev,jev,Ee) = J]^ V^dim(pe^/) dim{pev) X 

e 

(?;',e,D)cr 

This is normalized as 

{p'ev, j'ev: -e\pev, jev, '^e) = H ^ PevP'.Jje.jL H ' ^^^'> 

ev e 

Note that jev is in the dual to pev if (e, v) C F, and in p^v directly if {v, e) C F. 

For the case where there are no degeneracies d, d' the coefficients Ek^dd' simplify to Ek. 
Thus we can directly work with the basis of H invariant operators 

y dim(Kj 

labeled by k with the property 

-fc^-fc' = 4fc' • (38) 

We then obtain the states 



IPej;; jevi ^e) 



n 



V^dim(pe,;') dim(p 

e 



{v',e,v)cT 



v/dim(/ce 



X 



(39) 



X-Dp^^, {9ev')jev'^£v''^^P'^'"'^ ^e)i^^ime^^Pev, ^e) rrieiev^ Pi^viOev) ievjev ) 
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Which are normahzed as 



ev e 

To implement the G invariance at the vertices we can simply contract with inter- 
twiners rjy G Inv((8)p*^) where p*^ is the dual representation if {v,e) C F and the usual 
representation if (e,f) C F. Contracting all these we obtain the state 



-e) = n W v^dim(pe„/) dim(pe«) X 

{v',e,v)cT 

This is normalized as 

{p'ev,V'v,^'e\pev,Vv,^e) = '^P-P™ H^^'' '^''^ 11 '^^^^'^ • ^^^^ 

ev V e 

For the case without degeneracies this again simplifies to the states 

^ {v',e,v)cr 

Which are normalized as 

{p'evyV'v,K\pev,Vv,ke) = JJ H '^^^ 11 '^'^-'^^ • "^^^^ 

ev V e 

This should be contrasted to the basis of projected spin networks which is given in 
terms of H intertwiners t as 



I ^ev-; ^V-} Pe) = Yl '^^^^ev,... n X 

V e 

X I{pe,kev)ln^^j^^Dp^{gyev')jev,j^^:HPe,kev')j^^,,m^^, , (45) 

and is normalized as 

{Kv,Cp'e\kev,tv,Pe) = YlSk,^k',^Yl{t'Jtv)YlSp^p'^ . (46) 



e 
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3.3 Trimmed Complexes and Projected Spin Networks 



If the neighbourhood of the boundary of C is of the form Fx [0,1] we can make contact 
to the projected spin network space L^(G'''"'='/iJ'^"'). In projected spin networks the 
subgroup invariance is on the vertices of the boundary graph, rather than on the edges. 
If the neighbourhood of F is F x [0, 1], every boundary vertex has an associated internal 
edge V X [0,1], by "splitting" this associated internal edge we can move the subgroup 
invariance to the boundary vertices. To do so we need a square root of the E function, 
with the same subgroup covariance 

E{g) = I dg'F{g')F{g'-'g). (47) 

In terms of the basis coefficients this gives fj^fj^ = e^. This means that whenever we 
have a term of the form 

I dgdg'f{gg')F{g)F{g') 
with g,g' E G we can reparametrize with g = gg' and obtain 

j dgd~gf{~g)F{g)F{g~'~g) = j dgf{g)E{g), 

thus if the E function defines a projector we obtain E = F. 

While F inherits the symmetries of E, that is, F{hgh~^) = F{g), for all h E H, we 
generally have that F{g~^) ^ F{g). Therefore the amplitudes constructed from F will 
depend on the orientations of the faces, and we will need to keep explicit track of the 
orientation of the group elements associated to the edges in question. We do this by 
writing g^^v' = fi'^L for {vev') C / instead of gef, and we will split these as 

Consider now the partition function associated to a "cylindrical" 2-complex C = F x 
[0, 1] with boundary equal to two copies of F, F^ and F^. For simplicity for this section 
we will specialise to the case u = 6. It is straightforward but notationally cumbersome 
to extend the discussion to the general case by introducing a square root of the face 
weight. 

This partition function can then be factorized into the operators defined by 



v£e' v£e',e'£f (vev')cf 

X n ^i9ve9vev'9ev'9v'e'9le'9v'evgin^9e"v)E{g^ev')F{gl^,)F{g{„^){4S,) 
e&Ve 

where {e"v) iye), {ev'), {v'e') and {yev') C /, and e G Fg. e" and e' are the edges vx [0, 1], 
and v' X [0, 1]. The group elements in the delta function are those corresponding to the 
half of the face / near the boundary edge e. Note that the group elements g^'ev are 
reversely oriented with respect to the boundary, they are on the "opposite side" of the 
half face in the delta. 
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fx defines a map from tlie projected spin network space 
to tlie universal boundary space 



via 



(49) 



ser 

Tliis is constructed sucli tliat we liave 



eer 

(50) 

Tliis can be seen by explicit calculation, liowever, tliese calculations are greatly fa- 
ciliated by the graphical notation we will introduce in the next section, and we will 
illustrate them using examples there. 

In the spin network basis for the case without degeneracies the /x map can be expressed 

as 



ev V e3v e V \rej 

In general, if the boundary of the 2-complex is of the form Fx [0, 1] we can factorize 
its spin foam amplitude into the amplitude on the "trimmed complex" Ct, and fi for the 
boundary graph. The trimmed complex is the complex with "half of the boundary faces 
taken off" , or, more technically, where we consider the boundary edges and vertices not 
to be part of the edge set and vertex set of the 2-complex but to be in a seperate set of 
boundary edges. We thus have C^, C^, Cj, T^, Fg as separate spaces, however, still with 
adjacency relations and orientations amongst each other as before. 



(52) 



with 



2^{Ct)[9ef] 



Yl'^9ef 

\eCf 



\ vCe I 

\eeCe / 



( 



X 



n 



n 



, eeCe,C/ 



(53) 
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Note that this crucially depends on the orientations of the faces touching the bound- 
ary, and the composition of amplitudes only has the natural interpretation in terms of 
combining complexes if the orientations match up. 

We see that /i embedds the projected spin network state space into the, in some sense 
larger, universal boundary space defined above, and clearly 2^{C) lives in the image of 
fi^, thus we can equally well see the partition function as an element of the projected 
spin network space. 

3.3.1 Subgroup spin networks 

fi will generically have a non-zero kernel depending on E' or F respectively. Thus we 
can actually see Z as living in the coimage of n, which will allow us to go to subspaces 
of projected spin networks, for example H spin networks, as was discussed in [52]. This 
the case for example in the EPRL spin foam model. 

We can realize this restriction to subgroup spin networks explicitly if there is an u' 
with the properties of u such that the E function satisfies 

E{gMgg'g"g'")E{g") = [ dhuj'{g' gh)E{g)uj\g" g"'h-^)E{g"). (54) 

We can then glue via subgroup integrations on the boundary edges. Again, it is easy to 
see that this leads to the correct gluing using the graphical notation in the next section. 
We assume that every face has at most one boundary edge and replace u on those faces 
with a;', as well as reducing the boundary group element to live in the subgroup. The 
partition function then becomes 




where u* is u if the face does not contain a boundary edge, and u' if it does. 

3.4 A graphical notation 

It is illuminating to illustrate the structure of the convolutions in the partition function 
using a graphical notation. This will allow us to explicitly keep track of the way group 
elements in different faces are identified. The graphical notation will have three ingre- 
dients, corresponding to the face amplitude w, the insertion of E around a face, and the 
insertion and integration of the gev 
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E{9-') 


F{9-') 



Figure 1: Ingredients of the graphical notation. 
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We will indicate w by a solid line. White circles crossed by the line indicate group 
elements that are multiplied together to form the argument of u. The fins of the line 
indicate the order in which the inserted group elements should be multiplied. Two circles 
joined by a dotted line indicates the same group element. If one of them is crossed the 
group elements should be inverse of each other. The E function is indicated by a grey 
circle, the cross indicates that the argument in the E function should be the inverse of 
that it is linked to. We will represent its convolution square root F with a grey half 
circle. A white box indicates an insertion of oriented group elements into the lines that 
pass through it. This is given in Figured! 
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Figure 2: Some possible relations among the ingredients. The upper identity holds for 
u = 6, the second is a consequence of the definition of F, the third is the 
subgroup property. 

Using these ingredients we can represent the structure of a face containing four edges 
as in Figure [31 The structure of two half faces, and their composition in the projected 
spin network space, is given in Figure |H It is now a straightforward application of the 
relations illustrated in Figure |2] to see that the two half faces do indeed compose to a 
full face. 

We can write out the composition illustrated in HI Let that face be / = (012345670), 
with 0, 2, 4, 6 G C„ and 1, 3, 5, 7 G Cg. Taking care of the orientations in the delta function 
and the F, the integrand for a face on the right hand side reads 



'5(fi'0l5'0125'l2fi'235'23fi'210fl'70fl'70)5 (fl'455'4565'565'67fl'67^6545'34fl'34) X 
xE(^7012)F(^?2^3)^(^770)^(^?456)i^(^76^7)^(^?^4)'5fe0^674). 



(56) 



After combining the delta functions we obtain: 
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Figure 3: The structure of a face. 




Figure 4: Two half faces convoluted. 
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5 (fi'0lfi'012fi'l2fi'23fi'23fi'34fi'34fl'45fi'4565'56fl'67fl'67fi'70fi'70) X 

Which, with equation (1471) . gives the amphtude for a face in the normal partition 
function, pictorial represented in figure [3l 



^ (90190129129239234:934:94:594569569679670970) X 

xE{goi2)E{g234)E{g456)E{967o). 



(58) 



If we further have the subgroup property on E we can replace the convolution in 
the G with a convolution in H, thus reducing to the sub group spinnetworks. This is 
illustrated in Figure [51 



-o- 



I I 



I . I 



,--0-- 



-o--' 



I I I , I 



Figure 5: Two half faces convoluted using the subgroup property. 



4 Gluing of 2-complexes 

We will now show how the different Hilbertspaces discussed above give correspond to 
different gluing operations on the 2-complexes. To simplify illustrations we will focus on 
the case of the 2-complex dual to a triangulated surface, and various related 2-complexes. 
In that case the dual 2-complex has trivalent vertices and bivalent edges, see for example 
Figure [6l The 2-complex contains a central face with four edges, which has the structure 
illustrated in Figure |3l 

The partition function for this part of the amplitude is given in [71 This ampltiude 
can be obtained by various gluings from different building blocks, depending on how we 
take it apart. 
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Figure 6: A part of a triangulation (dashed lines), and its dual (solid lines). 




Figure 7: The part of the partition function corresponding to the part of the triangula- 
tion shown in Figure El 
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4.1 Face to face: the universal boundary space 

The first gluing is simply that corresponding the universal boundary space. As noted 
above, arbitrary 2-complexes can be glued along arbitrary edges. Thus we can in par- 
ticular simply take the faces of the 2-complex as individual partition functions. These 
partition functions for single faces were introduced in the companion paper as effective 
face weights Uf. 

Uf = Z{f) (59) 

We obtain one such Uf per type of face. In particular we have only one such effective 
face weight for a regular complex. These can then be composed simply by equation fl27p 
to yield arbitrary 2-complexes. This is sketched for the complex above in Figure [Hi 




Figure 8: The composition of the 2-complex from faces. Each small double arrow indi- 
cates a composition in the universal boundary space associated to one edge. 

Note that the composition in l-L\fBS given by (!27|) can be localize to a subset of the 
boundary graph. That is, given F G C and F' G C and a graph F** that is a subgraph of 
F and F', we can treat only the common sub graph as boundary, and glue amongst it, 
yielding, 

tr^' Z^{C)Z^\C') = Z^^'r^^'iC Ur, C), (60) 
where F Upb F' is the graph obtained by identifying F** G F, F' and then deleting the 
identified graph, and tr'"'' indicates integrating the elements of T-i\jBS ® '^'ubs against 
the element Yl^^S{gevg'ev~^)- 

Note that due to gauge invariance at the two-valent vertices Uf depends only on as 
many variables as / has edges. It will be convenient to also introduce u'j- such that 

^f{9ve, Qev', gv'e', de'v", ■ ■ ■) = 1^/ (fi'veS'ei)' , fi'D'e'5'e'D" , • • • )• 

4.1.1 A special case: Wedges to wedges 

In the case of the complex dual to a triangulation each type of face can ocurr, thus 
effective face weights are not a convenient choice for constructing theories. However, 
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we can construct a second kind of dual 2-complex C constructed from so called wedges. 
That is, we take as faces the intersections of the faces of the dual 2-complex and the 
simplices. In our 2-dimensional example this means that each triangle now contains 
three such wedges. These can then be composed in the universal boundary space again. 
The advantage is that now we only need one type of amplitude that we are gluing, the 
content of a simplex. This is illustrated in Figure [9l the new complex C is on the right. 




Figure 9: The composition of the 2-complex from the content of the simplex. Each 
small double arrow indicates a composition in the universal boundary space 
associated to one edge. 

Calling the complex of wedges a* , and using the trace from equation (160|) we can thus 
write the entire partition function for a dual complex made from wedges, 

Z{C') = iv^ ^ Z{a*) (61) 

V 

4.2 Half face to half face: the spin network spaces 

As noted above, the gluing of two trimmed partition functions in 'hJpsN ^^^o generates 
a natural composition, 

lz^iC),Z^iC')S' =Z(COrC'). (62) 

\ / PSN 

Due to the topological restrictions near the boundary required for trimming the parti- 
tion function, all edges of F in CUrC are bivalent, and CUrC is the 2-complex obtained 
by identifying the two copies of F and then erasing the bivalent edges in Fg and then 
the bivalent vertices in F^. This is the composition usually done in spin foam models. 

For the case of a 2-complex dual to a triangulation this is a natural type of gluing, 
an example is illustrated in Figure [TUJ The half faces and composition on the left hand 
side are exactly those illustrated in Figure HI 

4.3 Half faces around a face 

We can generalize the above gluing by using the partition function with 6 functions on 
the half faces, and gluing around a face amplitude u in the following sense. 
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10: The composition of the 2-complex from half faces. Each small double arrow 
indicates a composition in the projected spin network space associated to one 
edge. 



Figure 11: The composition of half wedges around a face. 
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For a set of boundary edges e^, a = 1 ... n in a set of trimmed complexes Ca, we 
can form the new partition function on a trimmed complex with an interior face / = 
{v,ei,v',e2, ...,en,v) 



This is illustrated in Figure [TTl 

In this way we can parametrize the partition function of the dual complex of a trian- 
gulation by the complex of trimmed wedges in a simplex and the face weight. 

5 Standard spin foam models in holonomy language. 

In this section we show that the BC, EPRL and FK model can be expressed in the 
language above. To do so we will relate the Holonomy model defined above to the 
operator models of [TB]. We can then give natural generalisations of the BC, EPRL and 
FK model to arbitrary, and, in particular, finite groups. 

5.1 From operators to holonomies 

We can relate the OSFM to the ones by insertion of the group integrations gef at each 
pair of edge and face by using the orthogonality of group elements: 



The structure of the manipulation is easiest to see using an extension of the graphical 
notation. We will decorate the lines with irreducible representations. We can then break 
them into segments with indicies living in the representation with arrows indicating the 
ingoing and outgoing indices. In Figure [12] we use this calculus to express equation 
( 164|) . A line joining two objects indicates the composition of tensors. A line coming into 
a tensor indicates an index in a representation space (downstairs) and a line outgoing 
indicates an index in the dual space (upstairs). 




(63) 




(64) 




b 



b' 



Figure 12: Orthogonality in graphical notation. 
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We will start by deriving the holonomy formulation given in the previous section, 
from the edge operator formalism of |18]. There we have, for a colouring of the faces 
by irreducible representations p, an operator P on each edge. These operators are then 
contracted according to the structure of the 2-complex. For simplicity we will focus 
on the case where the edge orientation and the face orientation agree. Then we have 
Pe £ Hom(p/ ® pf> . . ■ Pf") where /, /', /" 3 e. In all models studied so far except the 
KKL model [671 ESj P has an additional factorisation property which implies that it can 
be expressed in the following way: 

Pe^$) = PgU Ell ® eX 4c PgS2) , (65) 

where (a) is a multi index ranging over a/, with f 3 e, and Pq is the projector 
on the gauge invariant subspace. The E G L{pf) are a set of one linear operator per 
representaton space which satisfy relations ensuring that the fiducial orientations do not 
enter in contracting the edge operators. These are not group covariant. We will usually 
supress dependence of E on the representation space as it should be clear from context 
on which space it acts. If we write this projector explicitly as a gauge averaging we 
obtain the following: 

^e(a') = / dgvedgev ^py fo^^ (S't.e) ^ Ep/J^{gev)- (66) 

The graphical representation of the edge operator is given in Figure [121 



= E 

Mpi) I I I I ^ 

tTTt 

Figure 13: Factorisation of the edge operator in graphical notation. 




Inserting the resolution of the identity flMl) . and introducing a sum over representa- 
tions, we obtain 

Pe^,)= dg^edgev' Yldgef i^Dp^l^^^g^e) Dpj[f^{gef) Dp/J^{gev,)x 

f^e f3e 

Y.dha{p,j)DpJ^{g:})E!^^. (67) 



X 

Pef 
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with Yip j: ^^'^{Pef)Dp^jJ {g^j) Ej being the function E. The representation matrices 
around a face can be contracted to a character, and as we decoupled the representation 
label on the operators E the sum over pf can be performed exactly to arrive at the 
distribution ([3]) and we arrive at (|2]). Note that the conditions on E that imply indepen- 
dence of the fiducial orientations now imply E{g) = E{g~^) and the end result is indeed 
invariant under reversing the orientations. 

The graphical representation of this insertion of identities is given in Figure [HI 




Figure 14: The E functions. 



In the companion paper [15] as well as in [12] more general edge operators that do 
not have this factorisation property are also considered. Then there still is a holonomy 
formulation, given that the face amplitude is simply the dimension of p/. We briefly 
recall this construction and illustrate it in the graphical calculus in Figure [151 

The dimension as a face amplitude can be given by the trace of the identity operator, 
or a closed circle in graphical notation. The crucial ingredient is then a function on n 
copies of the group, given simply by 

. . ..gef) = (ll^MPef)DpJ;ig^j)\ P^y (68) 

The reason for this formula is clear by Figure [151 From this diagram we can also 
see that the C functions are glued with delta functions. In particular the integrand is 
schematically given by 

\{C{glfgtf\...)\{ n %eVeV')- (69) 

e V {eve')Cf 

In GFT language the C are simply the propagators and Y[(eve')cf ^^9ef^e'f~^) 
vertex function, or interaction, that glues them. 

5.2 The specific models for Spin(4) 

Having in hand the relationship to the operator formalism, we can now easily give the 
formulation of the various established spin foam models in the holonomy language. We 
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Figure 15: The C functions. 



begin with the BF, BC [23] and EPRL model [53] which are quite straightforward, and 
then discuss the FK and BO [2T1 [2D] in the subsequent subsections. 

5.2.1 BF, Barrett-Crane and Engle-Pereira-Rovelli-Livine 

We specify now to G = Spin(4) and H = SU(2) = Spin(4)djag the diagonal SU(2) 
subgroup, with irreps labelled by p and k respectively. The operators E in L{p) for the 
first set of models are given by the following: 

Ebf = 1, (70) 

for BF theory, 

Ebc = de(p) ^(P' 0)^(P' 0)^ (71) 

k 

for the Barrett Crane model, where de(p) is an arbitrary edge meassure factor, and for 
the EPRL model [Ml ESI SB] we have: 

Eeprl = de(p) J2 ^ (P' Hp, k)I{p, k)^ (72) 

k 

where we write p-yi^k) = (^^^k, ^-^-p^k^. 

We can then obtain the E{g) by the formulas of the previous section. For the EPRL 
model the functions E^^pj^^{g) is given by 

EIprM = dim(p)d,(p)5 (p, p,{k)) tr, {D,,{g) lip, k)I{p, k)^) , (73) 

p,k 
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and for Barrett-Crane we obtain 

Ebc{9) = Yl dim{p)d,{p)6ik)6 (p, {k', k')) tr, [D^{g) I{p, k)I{p, k)^) . (74) 

p,k,k' 

The SU(2) irrep k' plays a very different role than in the EPRL model, namely it 
restricts the form of the irrep p but does not appear in the injection maps /(p, k). We 
can further simplify this by noting that 

tr(fc,fc) {D^k,k){9) I{{k, k), 0)I{{k, k), 0)t) = I{{k, k), OYD^kMHik, k), 0) 

= tTk{g+g—')dim{k)-\ (75) 

with the Spin(4) element g decomposing into the left and right SU(2) as g = {g'^,g~). 
Thus for choice de(p) = 1 we have simply 

Escig) = Si9'^9~~^) = ^Spmi4U,,{g), (76) 
where, for a general subgroup H C G we write 

Snig) = [ dhSigh-') (77) 
Jh 

for the delta function that force a group element to lie in the subgroup. 

Thus we arrive at a particularly simple form for the Barrett-Crane model as an integral 
over a product of SU(2) delta functions: 

2bc{C) = I {U^h%^ ij^-f^^") (^^ (^")'')) ill^^3j)S{gj?j . (78) 

For the FK model the E function is most easily expressed in terms of coherent states, 
that is, the eigenstates of the Lie algebra generators. For an SU(2) representation 
labelled by the half integer k these are the states a{n) that satisfy (n-L)afc(n) = ikakin), 
where n is a unit vector in M^. All afc(n) for the same n differ at most by a phase. The 
E function for the FK model is given by 

Efk{9^.9-) = X;dim(A:) f dn (^(n)/}, (<?+)«, (n)) (4(n)Dfc((7-)a,(n)) , (79) 
k ^ 

which is well defined as the phases of Q!fc(n) and al{n.) cancel. 
BF theory of course is simply given by setting 

EBFi9) = Sig). (80) 
The coefficients in the basis of section 12.31 for the various models then are: 
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• BF: el = 1 

. BC: e^ = de(p)5MEfc'5(P,(^',^')) 

. EPRL: = d,(p)5 ((i±lfc,^A:) ,p) 

We see from the conditions derived in section [2731 that for real de(p) all these spin 
foam models are indeed real holonomy spin foams. 

5.2.2 Freidel-Krasnov for 7 > 1 

The FK model is defined in terms of coherent states. For 7 < 1 it is equivalent to the 
EPRL model. To obtain its basis coefficients for 7 > 1 we need to work some more. In 
this section we use the shorthand \m) = a{m). 
The operator Ep for the FK model is given by 

E„ = de(p) /— |m)2^'^(m|2j'' ® |m)2j''(m|2^'~ (81) 
J 47r 

with the representations p = subject to constraints similar to those in the EPRL 

model. Let 

i?,= |(^fc,l^^) :^ = 0,i ...}. (82) 

We can write 

E{g) = 5^dim(p)d,(p)5«^(p) f ^{m\g^\m)'^^ {m\g+\mfr . (83) 



Now we want to compute 



,-^^TTl\p,k)E,I{p,k) (84) 



Note that /(p, k) is given by its matrix elements in the coherent state basis, thus if 
p = with = -^k and j"*", j~, k admissible we have that 

(/(p, k)\ny',\n+f^+ ^\n-f^-) 

= "^^^^ \ r2,n+)^(^^+^-^-^-')(n,n")^(^+^-""^-")e(n+,n")^(^--+^-^-^) ^^^^ 
with e(-, ■) being the invariant bilinear form and 

^2 ^ (j+ + r + ky.jj- + j+ -k)\{k+ j+ - j-y.jk + r - ry. 

(2A:)!(2j+)!(2j-)! ^ ^ 

is derived in the book by Kauffman and Lins 
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Since /^(p, k)EpI{p, k) is proportional to identity, we only need to find 

{m'\Pip,k)E,Iip,k)W2f>'), (87) 
for the case p E R^. Then this is equal to 

4^i±l f ^(i/2,m)2('=+i+-r)(i/2,,ri)2('=+^--^-+)e(m,m)2(.^-+i+-^) 

^ J 47r (88) 
(m, V2)2(^-+^'+-^'")(m, i/2)2('^+^'"-^^+)e(m, m)2(^'"+^'+-'=) 

But we know ^ ^ 

|(V2,m)p= (^cos0 , |(V2,m)p= (^sinQ (89) 

where 6 is the angle between the direction of m and the north pole. We can introduce 
polar coordinates 



e 



^ = ((y2pMt(p,fc)4j(p,fc)|V2)^^ 



2k 4- 1 W2 _ _ ^^^^ 

= 4 — ^27r2 / c/r/2sinr/cosr7(cosr/)2('=+^'^-^'")(sinr7)2('^+^'"-^'^) 
4vrC^ Jo 

^ 2fc + 1 (fc + j+ - j-)!(fc + J- - j+)! 
" C2 (2A; + 1)! 

We thus obtained 

' 'u*+]- + m3-+i*-k)\ ^ ' 

Substituting = ^^k we have that for p = the basis coefficient is given by 

= d ((7 + l)fc)!((7-l)fc)! .92) 
' {jk + k)\{jk - k)\ 

For 7 > 1, is the minimal representation in the decomposition and 

el = 4 (93) 

Finally, making the conditions on p explicit again we can write 

Thus the FK model also falls into the class of real holonomy spin foam models. 
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5.2.3 Baratin-Oriti 

Let /3 = 2^. We will use the following 7 dependent transformation from SU(2) to 
SU(2): 

SU (2) 3 cos 6 + ina sin 6 = u = cos 6[s + m^^a sin (94) 

where the class angle 613 and the unit vector are determined by conditions [21] 

sin^^ = sin0, sign(sincos6'/3) = sign(cos6'), n/3 = sign(/3)n . (95) 

Let us notice that 

(n/^)-i = (u-y, (gug-y = gu^g-' . (96) 

The fusion coefficients for the BO model are given in [21] (eq. 37,58). From these one 
can derive 

- = du (u~^)D-^^ ((u^)~^) 

{m+in~)(m+m~) j m^n^^ ' m+n+VV / / 

I du Dil,^4u)DC,^4u') (97) 

Thus the E function 

Eig\g-) = di^O'^) dim(j-)^;;i- -)(^.^-)/^£.^.(^7+"^)/^l-™-(^?-"') (98) 
i+i" 

is equal to 

E{g+,g~) = dim(j+)dim(j-) J dudu Xri.u"^^9~'^)xA{^^y^^^ 9^'') (99) 

dudu 5{u-^ug-~y{{u^)-^u^g+~^) (100) 
We can simplify the equation by solving the first delta function for u so that u = ug^ 

E{g^,g-) = fdu S{{u^)-\ug-)^g+-') (101) 



In the case of /3 = 1 this i?-function reduces to the one of the BC model. The BO model 
is also a real holonomy spin foam model. 

6 The transfer operator for holonomy spin foams 

So far we have discussed at length the structural aspects of holonomy spin foam models, 
their boundary Hilbert spaces and their gluings. This raises the question if, having the 
boundary Hilbert spaces at hand, we can define a Hamiltonian dynamics refiecting the 
one defined by the spin foam models. 
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The first step in deriving such a Hamiltonian dynamics from a given partition function 
is to obtain transfer operators. In standard lattice theories, for example lattice gauge 
Yang Mills theory, such transfer operators correspond to finite time steps. To obtain 
the Hamiltonians as infinitesimal time evolution generators for such systems one would 
have to take the limit of infinitesimal time by scaling the coupling constants in time - 
and space directions in a certain way defined by the dynamics of the system [691 EO] 

The issue of obtaining the Hamiltonians is more involved in gravitational systems, 
as the lattice constants and the time separation are rather encoded in the boundary 
states of the system. Furthermore the question of exact and broken diffeomorphism 
symmetry comes in [T3lllllll6]. Only in the case that exact diffeomorphism symmetry is 
preserved in the discretization, can we expect the appearance of Hamiltonian constraints 
in the canonical formulation [HI [13]. If this holds also for the partition function the 
transfer operator is a product of projection operators from which the Hamiltonian and 
diffeomorphism constraints can be read off [771 ESI [IS [IS]- In this case no limiting 
procedure is necessary. If diffeomorphism constraints are broken, one can either attempt 
a limiting procedure involving the boundary states or alternatively attempt to obtain 
an improved model by coarse graining which then carries a notion of diffeomorphism 
symmetry [131 [II iH iSl SSI [13 113] • 

A third possibility is to adopt the view point that the dynamics is inherently discrete, a 
viewpoint which is for instance emphasized in the framework of consistent discretizations 
[581 [59] . In this case the transfer operator can only defined for finite time steps, and a 
limit cannot be taken (in general). 

Here we will consider the finite time transfer operator and comment more on the issues 
of taking the limit to obtain the time evolution generators afterwards. We will consider 
a space time lattice with a regular slicing in time direction, i.e. each (thick) time slice is 
of the form Tg x [0, 1]. The discussion can be generalized to some extend to an irregular 
lattice and a notion of local time evolution, see for instance [HI [19] for a discussion in 
classical Regge calculus. 

The definition of the transfer operator requires a choice of slicing of the underlying 
lattice and initially we choose one which will make the transfer operator as similar as 
possible to the one encountered in lattice gauge theory [691 IZOl [951 [IS]- As we will see 
such a slicing fits well to the universal boundary Hilbert space introduced in section 13.11 
In this formulation the effective face weights introduced in [I5] and recalled in section 
14. II will play a prominent role. From a simplicial geometry view point the faces are dual 
to the bones of the triangulation, which carry the curvature. Thus this slicing offers a 
new perspective on the Hamiltonian dynamics and the semiclassical limit: it does not 
concentrate on the vertex (i.e. simplex) amplitude but on the gluing of simplices around 
the bones, where the curvature and hence the essential dynamical information resides. 

On the other hand this slicing is somewhat unusual in discrete gravity, where one 
often builds a transfer operator by gluing simplices to the hypersurface [H [IHl [IH]. In 
this case equal time hypersurfaces can be understood as (dual to) {D — 1) dimensional 
triangulations. As we will see in the course of the discussion we can switch to a slicing 
more adapted to a simplicial viewpoint by using the fi map between the projected spin 
network and universal Hilbert space introduces in section 13.31 The /i map can then be 
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understood to project onto the solutions of the simphcity constraints - the (stripped) 
transfer operator will be sandwiched between such projectors. 

6.1 Transfer operator for general models 

Transfer operators can be defined if we have a discrete "time" direction in our 2-complex 
in the sense that the complex C is the 2-skeleton of the complex [Cg x [0, 1])", with 
C" X [1] = C^"*"^ X [0]. In that case we call the edges and faces in the various C" spatial, 
and the other edges and faces temporal. 

We call the graph of horizontal edges F^. Then the partition functions Z{Cs) and 
Z{V^ X [0, 1]) act naturally on the space ^^BS' ^^^"^ simply by multiplication. 

It follows directly from the gluing of in ^^BS ^'^'^'^^'^i^S (127|) that the partition 
function of C can then be written as 

Z{C) = {Z{Cs)Z{Vs X [0, 1]))" Z{C,) (102) 

This has the structure of a partition function written in terms of transfer operators. 
For the rest of this section we call Z(C,) = Z(F, x [0, 1]) = K and Z(C,)Z(F, x 
[0,1]) = T, that is 

T = WK, (103) 

and 

Z(C) = T^'W. (104) 

The operators W and K can be written very efficiently in terms of the effective face 
weights Uf which we recalled in section H?T] in fl59l) . Recall that these are simply given 
by the amplitude of a face, as illustrated in Figure [31 

As the entire 1-skeleton of Cg is in the boundary space the only integrations in W 
are those involving the Qef. These are exactly the integrations one performs to obtain 
the effective face weights Uf. The operator W acts as a multiplication operator in the 
holonomy basis of the universal Hilbert space and is just given as a product of the 
effective face weights 

w[9e.] = n ^AiaV)^ (105) 

feCs 

where the notation {g}^ indicates the set of all group elements belonging to the face /. 

For K we have a similar simplification. Its structure is illustrated in figure [161 In 
addition to the integrations over the Qef we have also to take into account the integrations 
over the g^^ associated to the time like edges, which we will denote by gl^. By g'l^ and 
(7g^ we denote holonomy variables associated to edges in the spatial hypersurfaces at 
two consecutive time steps. The integration kernel of the operator K in the holonomy 
representation is given by 
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(106) 



where each effective face weight u!'j:{{g}^) depends on four variables gl^ associated to the 
two time hke edges of the time hke face, a set of variables g^^ associated to the space 
like edges of this face shared with the graph Tg x [0] and a corresponding set g'l^ shared 
with X [1]. 




Figure 16: The composition of two effective face weights as it occurs in the definition of 
K. 

Each time like edge has one boundary vertex in x [0] and one in F^ x [1], and each 
vertex in F^ x [0] or F^ x [1] has only one temporal edge going out. Thus we can drop 
the e in gl^ and simply write g^. 

Recall that the effective face weights Uf naturally live in the universal boundary space 
of a segmented hne, and thus have a gauge freedom that acts as gye — > gvdve- Hence we 
can apply gauge transformations to the vertices in F^ x [0] and F^ x [1] such that all the 
group variables associated to the time like edges are equal to the identity. This allows 
us to write the integration over these group elements associated to the time like edges 
as gauge projectors on the vertices of the graph F^ 

K = PgKoPg. (107) 

where 

PcW^gL] = / ( n ^U^^^'evdidD), (108) 

is simply the projection on gauge invariant functions on the vertices of the graph. 

Ko is obtained by setting the group variables in cu/ or equivalently u'j associated to 
the time like edges equal to the identity element 

Ko[g'l^,gU = Yl^f(e){9v,e, gev2, 9v2e^ g'evi) (109) 

e 

where Vi,V2 are the source and target vertex of the edge e and /(e) = e x [0, 1] is the 
time like face associated to the edge e G F^. 

This structure is illustrated in figure [171 A simple application of the top relation in 
figure [2] will return us to figure [161 



32 




Figure 17: The decomposition of the effective face weights into reduced face weights and 
gauge projectors. 

The operator Kq does not map from the universal boundary space to the universal 
boundary space as it does not produce a state with G invariance at the vertices but 
acts in L^(Gl'"'="l) instead. The projector Pg then brings us back to the gauge invariant 
universal boundary space. 

It will be convenient to write Kq = Yle as an operator in the usual sense, in 
terms of left and right shift operators {L{g) > /)(■) = f{g~^-) and {R{g) > /)(■) = f{-g) 
respectively. To this end we introduce a bra-ket notation such that we write an element 

1(91,92) e L'^iG X G) as f{9i,92) = {9i,92\f)- 
The operator Ke can then be written as 

Ke = dgAd9BYl{dgiE{gi)) ujf{gB92^) Li{94)Ri{9i9a)L2{9b93)R2{92)- (^O) 
^ i=i 

where the left and right shift operators act as 

Li{ga)Ri{9b)L2{gc)R2{9d)\9i,92) = \9a^9i9b, 9c^929d) ■ (111) 
In summary the transfer operator is given by 

T=w-K= (n^/j iu^oj (n^^) (n^^.O 

where W acts as a multiplication operator in the holonomy basis and factorizes over the 
faces, Pg,v is the projector on the G-gauge invariant subspace at the vertex v and we 
defined the action of Ke in equation flllOp . 

6.2 Transfer operator in the spin network basis and simplicity 
constraints 

The operator Ke is simplest in the spin basis. As Ke just acts on one edge it is sufficient 
to consider the one-edge Hilbert space L'^{G x G). A basis for this space would be given 
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by |pi, p2, ^2,^2), where ^2, J2 are magnetic indices in the representations pi, P2 
respectively. 

However, to compactify notation it is useful to introduce a basis adapted to the H 
group in a given p representation. The basis we introduce is of Gelfand-Tsetlin type 
[62]. It is labelled by 

j^{k,d,m}, \{k,d,m}) = I{p,ky\m) (113) 

where k is the label of the H representation, d is the multiplicity index (in the multiplicity 
free case it will be omitted), and |m) is a basis in the k representation. Thus we 
replace the (four) magnetic indices of the p-representations in L^(G x G) by four indices 

ji ^ {ki, di, TTlj}. 

The basis in the space of H invariant functions is thus spanned by 

|pi, {ki, nil}; k; p2, {k2, ms}) = ^ V |pi, {ki, mi}, {k, m}; pa, {k2, ms}, {k, m}) 

Vdim(/c) ^ 

(114) 

the same functions as introduced in (l39i) . 

As follows from (11 101) and picture [T7] the operator Ji'eCan be decomposed into the 
following more elementary components 

= KIK'K^ KiKl , (115) 

where 

= j dgAdgs ^f{gBgA^)Ri{gA)L2{gB) = j dgujf{g)Ri{g) j dg' Ri{g')L2{g') (116) 

and 

KI = j dgi E{gi)Ri{gi) = j dg2 E{g2)R2{g2) (117) 

= j dg:, E{gs)L2{g3) Kl = j dg^ E{g,)Li{g,) . (118) 

Let us notice that K2 and Kl commute with the rest of the operators and KfK'^K^ 
commutes with Pg- 
We also have 

K' = ^uPPp (119) 
p 

where Pp is the projection onto the subspace in L'^{G x G) spanned by the following 
orthonormal basis labelled by ji, ja 

^===y' |p,ji,i;p,j,i2) • (120) 

Vdim(p) 
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The operators Kn^ act on the basis ^i, Ji; P25 ^2, ^2) as follows. The operator Kl 
changes only the index ji = {k, d, m} by multiplication of the matrix 



^p, dd^^^hk ■ (121) 



Similarly the operators fCj? -^35 -^4 act only on the indices j2,i2,H respectively. 

In the multiplicity free case we have a straightforward eigenfunction expansion of both 
K2KI and K^K^K^. The eigenvectors with non-vanishing eigenvalues for K^K^^K^ are 
given by 

|p, {h, mi}; {/C2, ^2}) = , ^ ej|p, {/ci, mi}, {k, m}; p, {k, m}, {/cg, "^2}) 

YEfedim(/i;)(e^)2 

(122) 

with corresponding eigenvalues 
iqK'KI\p,{ki,mi}-{k2,m2}) = dim(/i;)(ej)2 j |p, {/ci, mj; {/i;2, m2}) . 

These are also eigenvectors for K\K\ with eigenvalues 

. (124) 

In summary the eigenvectors with a priori non-vanishing eigenvalues for K^, are given 
by f ll22p with eigenvalues 

^--^AY^^^mf) • (125) 



(123) 



Thus the eigenvalues are independent of the labels mi, m2 and come with a multiplicity 
dim(/ci) dim(/c2)- Note that K^, vanishes on states "^i, ji; ^2, ^2) with pi 7^ p2 as well 
as on states with p = pi = P2 but orthogonal to (11221) . 



In the gravitational models the E functions and therefore the impose the simplicity 
constraints. Hence we can say the same of K^. - it maps onto a subspace of the universal 
boundary Hilbert space on which the (primary) simplicity constraints hold in some 
form0. As the map to a subspace which can be interpreted as solutions to the 
primary simplicity constraints let us also consider the question whether W , or some 
suitable subset of holonomy operators, leaves this subspace of the universal boundary 
Hilbert space invariant. 



"^As the discrete form of the primary simpUcity constraints do not even commute weakly one has the 
choice to impose them strongly, i.e. as operator equations, as in the BC model [24j or in a certain 
weak form as in the EPRL model [55] . 
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is a multiplication operator that factorizes over the spatial plaquettes. The contri- 
bution from a given plaquette is of the form 

(126) 

where the edge e joins v with v'. From now on we will consider only the multiplicity 
free case, thus omit the index d. 

These are contractions between Sp-^pSp^pe^, and basic holonomy operators 

V^pi.fci.mi ;fc;p2,fc2,m2 = -^p ifc.ml (^ef i } ko . .mo\(9 ev2) ■ (127) 

■m 

Note that ippj^,ki,mi;k;p2,k2,m2 actiug on the constant function creates the states 

|p, {ki,mi}; k;p2, {k2,m2}). 



Thus the action of W involves the multiplication of holonomy operators of the form 
f ll27p . We will therefore consider the product of two such holonomy operators. This is 
a straightforward calculation, in which one first rewrites 

D; ,{g) D^;, ^^{g) = C^'''"/C^^'']^:D^,^„{g) (128) 

p",i",j" 

with C-yy, the Clebsch Gordan coefficients of G. Here we assume that G is multiplicity 
free, i.e. there is maximally one copy of a given irrep in the tensor product of two irreps. 
We will assume the same property to hold for the subgroup if, furthermore we already 
assumed that there is maximally one copy of a given H representation/;; in a given G 
representation. 

In this case the Clebsch Gordan coefficient for G contracted with the maps /(p, /c) 
reduce to the Clebsch Gordan coefficients of H which we can write in the Gelfand- 
Tsetlin-like basis as 

f^pp'p" {k".m"} _ ^kk'k" m" (^0Q\ 

^ {fc,m},{fc',»n'} ~ mm' • l-^^^J 

Finally the summation over the index m in the holonomy operators (11271) leads to the 
following contraction of Clebsch Gordan coefficients 

Cmrn'm''^mm'm'' = ^(^5 k ,k ) 5 m" ,m" (130) 

m,m' 

where d{k^ k', k") = 1 if k, k', k" couple to the trivial representation and is vanishing 
otherwise. 

The product of two holonomy operators of the form (I127p is therefore given by 

^ p\,k\,m\ -jk; p2,k2,fn2 ^ ^Pj,fcj,m'j^ ;fc' ; P2i'^2''"2 

V c^'^^f ' „ eik, k', k") c^^fc^^ „ ^ „ ^„ ^„ . „ ^„ ^„ . (131) 

p'{,p'^,k",k'^,k'^ 
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Here we sum over all repeated magnetic indices. The result is again a linear combination 
of basis states f ll27p . 

From this expression (11311) we notice the following: (a) Even if initially the represen- 
tations satisfy pi = p2 and p'^ = p'^ this will in general not hold for the basis states 
appearing on the right hand side of fll3ip . (b) Consider the case that we multiply two 
basic holonomies (11271) satisfying p = p\ = p2 and p' = p\ = p'2 which have been con- 
tracted with and e^, in the k and k' index respectively. The basis states in the product 
holonomy are then contracted with 

Y^eiei B{k,}^ ,k") (132) 

in the k" index. In general the product of holonomy operators of the form (11271) will 
not generate a proper subspace. An exception is the Barrett Crane model, in which the 
^'-function has an enhanced symmetry that allows a restrictions to spin network states 
with = 0. In this case ^(0, 0, /c") 7^ indeed leads to the condition k" = 0. 

The fact that the holonomy operators (I127p do not lead in general to a proper subspace 
might not be a surprise to the expert as the secondary simplicity constraints, which are 
conditions on the holonomies |100[ [5|[5T | [50] are usually not imposed in spin foam models. 
Indeed the hope is that the imposition of the primary simplicity constraints on two con- 
secutive time slices leads to the automatic imposition of the secondary constraints, whose 
function in a canonical formulation is to ensure that the primary simplicity constraints 
are preserved under time evolution. Later on we will redefine the transfer operator to 

make this notion more explicit. Here we just note that one possibility is to consider 

1 1 

T' = Kq PqW PgKq as we can interpret Kq to impose the simplicity constraints. 



The operator W naturally factorizes over plaquettes and we can expect that it leads 
to the curvature term F in the gravitational Hamiltonian constraints, which are of the 
form FEE with E representing flux (infinitesimal shift) operators. On the other hand 
we can also seek an expression which factorizes over the vertices of Tg. Such a form 
brings as back to the usual vertex amplitude representation of spin foams and shows the 
consistency of the procedure. 

For this calculation the gauge invariance at the vertices of is essential, hence we 
choose the basic states \pev,Vv,ke) introduced in fH51) which, using the GT basis, are 
defined by 

U(ev) Vdim(p)e„ 



Pel)) 5 ^e) 



n.v/dta^ 

(133) 



Let us note that W preserves the gauge invariant TiuBS (with respect to the gauge action 
at the vertices of the underlying graph). 
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We want to compute the matrix elements 



K) . (134) 

ly is a multiplication operator in the holonomy basis, hence to compute the matrix 
elements we introduce two resolutions of unity into the matrix elements fll34p . These 
resolutions of unity lead to an integration over the group elements g^^e and gev2- 

The 

holonomy operator associated to a given edge e is then given by 

(W^e){i;,i^(,+i,}({P/}) = n X] Dpj{g^^e)i,.^,{kf,,mf,}el^^.Ppf{gev2){kf,,rnf,},ii(,+^ (135) 

where the contribution {We) come with magnetic indices {ifeif{e+i)\ that are contracted 
between the different edges of a face. (Relative orientation of the edges and the face is 
unimportant if the model is real, otherwise we assume here that these orientations agree.) 
The We are contracted and then summed over the pf (multiplied with CjPi dim(p/)) to 
obtain the full operator W . 

The computation then proceeds in the following steps which are completely analogous 
to the construction of vertex amplitudes in spin foam models, see for instance [TBI [M] . 

(a) The integration over the group elements g^^e and gev2 leads to the Haar projector 

Piiji-injniPu ■■■,Pn) = J dgDp^{g)i^j^ ■ ■ ■ Dp^{g)i^,j^ (136) 

on each half edge. For instance for the first half edge the projector is on the invariant 
subspace in the tensor product 

(b) The Haar projectors on each half edge can be split into a sum over a basis of or- 
thonormal invariant vectors or intertwiners rj of the corresponding representation space. 

Piljl-injniPU ...,Pn) = Yl \v)il-ir.{v\jl-U ' (138) 

V 

(c) The magnetic indices of the invariant vectors on the left in (11381) associated to 
half edges Vie and of the invariant vectors on the right in (11381) associated to half 
edges ev2 contract now at the vertices of the graph Tg. Thus all the magnetic indices 

of (We) 

{ifeif{e+i)} contract among each other. The magnetic indices ivie,jev2^ivie-ifev2 
associated to the representation matrices of our spin networks are contracted with the 
intertwiners riy,'!]'^ of these spin networks. 

(d) The resulting amplitude is the vertex amplitude A^^ for G-BF theory. Here a 
vertex of F^ is to be understood as a vertex in the following 2-complex: The edges e of 
Tg are 'horizontal' edges in this two-complex and labelled with interwiners rj^e, which 
appeared in the expansion of the Haar projectors. Additionally we have the spatial faces, 
on which W is defined and which are labelled by pf. There are additional 'vertical' edges 
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and faces, which carry the labels of the spin networks between which we compute the 
matrix elements. For each w G we have an edge pointing down and labelled with 77^ 
and an edge pointing up labelled by rjy. There are also two vertical (half) faces attached 
to each (half) edge {ve) e F^ which are labelled by p^e and p^g. The orientation of these 
faces is such that these agree with the orientation of e for the 'up' faces and are opposite 
with respect to the orientation of e for the 'down' faces. This 2-complex around a vertex 
is depicted in figure [T8] 



riv 




Figure 18: The labelling of the faces and edges around a spatial vertex. 

Thus A^^ depends on all the algebraic data attached to (half) edges and faces adjacent 
to the vertex v in this 2-complex. This includes the intertwiners ?7„e, flv-, fl'v which involve 
the representation spaces pf for faces sharing / as well as p^,e, p^g for edges in F^ sharing 
V. The vertex BF amplitude is then defined by the contraction of the invariant vectors 
in the corresponding representation spaces 

A^^iVev, Vv, V'v) = tfpe. trp'^^ tipf^p^,^,,, riyr]^> JJ r]^y. (139) 

ev 

(e) We are left with half of the invariant vectors arising from the Haar projectors associ- 
ated to the edges of F^. These are contracted with the /(p/, kfe)jf^mf^e'kfj\Pf^ ^fe)mfj'^^ 
part of the operators We as well as with the Y,m, HPvie, ke)i^^emjKPev2^K)m,i,^2 part of 
the spin network state | Pev,Vv, ^e) and the corresponding primed counterpart. Thus to 
each edge in F^ we associate the amplitude 

f p \{Pf}f3e,Pvie,p'y^e^Pev2,p'ev2 / \_ 
^ '^Hkef}f3e,ke,K V hie, '/eiij j — 

{'nvie\y{pev2,K)I{pvie, K)^ ® Kp'v^e^ K)I {p'ev2^ K)'^ ^ ^Kf^^Pf' kef)I{pf^ ^e/)^ ) \Vev2) ■ 

(140) 

(f ) Finally the matrix elements of W are given by 

{p'ev^V'v^KW\pev,Vv,K) = 
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where M collects all the dimension factors and face amplitudes 



Thus the matrix elements of W provide us almost with the full spin foam amplitude 
for an equal time slice. What is missing to obtain the full amplitude are the insertions of 

for the edge face pairs consisting of time like faces and space like edges as well as the 
edge amplitudes Pg for the time like edges. These are indeed provided by the operators 



6.3 The transfer operator on the projected spin network Hilbert 
space 

So far we discussed the transfer operator in the unsymmetric form T = W ■ K. In lattice 
gauge theory one often chooses rather a symmetric form Tlg = W^^'^KW^^'^. In the 
case of spin foams W might not be a positive operator. Indeed for the gravitational 
spin foams it is rather easy to construct a square root of so that we can consider 
T' = Yle^e^^PcWPcYlKy^ . We have seen that K^. is almost a projection operator 
and have given the eigenvectors and eigenvalues fll22lll25p . Indeed for the EPRL model 
the non-null eigenvectors (11221) reduce to 

P = (^^' ^T^^) ' "'i^' ^2}) (143) 

so that the index structure is the same as for (non-gauge invariant) H spin network 
functions, i.e. the spin network basis of the standard loop quantum gravity Hilbert 
space. By formally identifying the eigenvectors with this spin network basis we can 
define the transfer operator on the LQG Hilbert space. In principle this applies also to 
the Barrett Crane model, the non-null eigenvectors are however just labelled by one H 
representation label 

^;f^:=|p = (fc,fc),{0,0},{0,0}) (144) 

whose functionality is in addition quite different from the EPRL model. 

A more elegant and geometric method is to use the fi^s niap introduced in section [231 
which maps from the Hilbert space of projected spin network functions to the universal 
boundary Hilbert space. (As in this section we will therefore assume that the face 
weights are given by delta functions 00 f = 6g, however as mentioned there, this can be 
generalized.) The /i map was defined such that the effect of the time like plaquettes, i.e. 
K, can be written as 

K = firA.. (145) 
In the following we will suppress the index Tg. 



EPRL ._ 
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Thus we can define the transfer operator in the projected spin network space by 

Z{Ct) = {fi^Wfiy i.e. TpsN = fi^Wfi . (146) 
Remember that the integration kernel of fi is given by 

W ^{9vie9viev29ev29v2t9t2t9v2evi9tvi9tvi)E{9y^ev2)^{9v2t)^i9tvi) (147) 

where F satisfies J E{g)f{g)dg = J F[g)F[g~^g')f{g')dgdg'. F can be expanded in the 
same way as E and for the coefficients we have = /^/^. In the following we will 
assume that is a projector, i.e. = e^. Here t indicates a time like edge, i.e. g^t 
is a group element associated to a time like edge starting at the vertex v, gt^t a group 
element along a time like edge ending in v' . 

In the same way as for the operator K we can extract a gauge invariant projector 
from /i, i.e. 



/i 



PG\{^^e (148) 



e 



where /ie is given by 

{9vie, 9ev2\l^e\9viev2) = jd9 d9v^td9v2t 

S{9 vie9vi ev2 9ev2 9v2t9v2evi 9tvi )E{g..e.2)F{9:2t)mj ■ (149) 

In a spin network basis of the universal boundary Hilbert space and the projected spin 
network space respectively we obtain 

(pi, {ki, nil} ■,k; p2, {A;2, "i2}|Aie|p', {k[^ m[}, {k'^, m'J) = 

^^|^^Pip"^P2P' 4 fkMk[ fC^k2k', ■ (150) 

where we also used a Gelfand-Tsetlin basis for the non-gauge invariant projected spin 
networks. 

Thus the image of /ig is spanned by the non-null eigenvectors for Kg discussed pre- 
viously and given in fll43p and fll44p for the EPRL and BC model respectively. The 
co-kernel of /ie or image of in the projected spin network space is labelled by the 
same indices, i.e. 

{)f^ = |(A;,fc),{0,0},{0,0}) 
%™2= (^-^kj-^^k) Ak.mi},{k,m,}) . (151) 
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Thus we can for instance for the EPRL model formally understand the transfer oper- 
ator as an operator on the LQG Hilbert space. This might however not be very useful, if 
one wants to understand the structure of the transfer operator in terms of holonomy and 
flux operators. In particular the holonomy operators appearing in the transfer operator 
are G holonomies and thus act on either the universal boundary Hilbert space or the 
projected spin network space, both of which are Hilbert space over copies of the group 
G. 

6.4 Example: The BF model 

Let us first consider the BF model, see also [16] for a discussion of the corresponding 
transfer operator in the context of standard lattice gauge theory. BF theory will be the 
only case where the transfer operator will be actually a projector. 

In the case of the BF model we have = 9{p, k), where 6{p, k) = 1 ii k appears in 
the reduction of p over H and 9{p, k) = otherwise. Also, as the face weights are given 
by delta functions on the group we have = 1. 

Ke projects onto states 

'^p,fci,mi;fc2,m2 

^ ^/dim{k)9{p, k) |p, {/ci, mi} ; /c ; p, {/ca, ms}) 

k 

= lHp,ki)m^hSj^iJ{p,k2)j2m2\p,h,ji;P,i2,j2) (152) 

ii,ji,i2,j2,m 

which are just linear combinations of spin network states that have been subdivided 
with a trivial two-valent edge. The eigenvalues of K^. are given by 

Ap,fci,mi;fc2,m2 = ^^T^^^ ^(P, ^i) ^(P, ^2) ^dim(A;)6'(p, /c)^ = 9{p, ki) 9{p, k2) , (153) 

i.e. are equal to one as long as the indices in fp,fci,mi;fc2,m2 define a non- vanishing vector. 
Thus Kf> is a proper projection operator and acts as identity on states that can be 
embedded into the standard lattice gauge theory Hilbert space for G. 

The operator W is a multiplication operator in the holonomy basis given by the 
effective face weights associated to the spatial plaquettes. As E{g) = 6{g) the effective 
face weights are also given by delta functions evaluated on the holonomy around the 
plaquette 

■■) = % viei9eiV29v2e2 ' ' ' ) 

(154) 

Hence is a projection^ onto the states satisfying the flatness conditions. 

Hence the transfer operator T is a projector implementing the Gauss constraints, 
which impose gauge invariance, as well as the flatness constraints on the plaquettes. 

^As it involves the delta function it is not a proper projection operator in the case of Lie groups. A 
mathematical clean description can be obtained by interpreting W as rigging map which maps to 
the dual of some dense subspace of the Hilbert space, see for instance [77] . 
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6.5 Example: The BC model 



Let us also discuss the Barrett Crane model as here an enhanced symmetry of the E 
function allows to reduce the boundary Hilbert spaces. 

For the Barrett Crane model the coefficients factorize = SkflJ2k'^p,(k',k')- The 
operator Ke projects onto states 

Vk' = \{k',k'),{0,Oy,0;ik',k'),{0,0}) (155) 

and the corresponding eigenvalue is given by A^/ = (dim(fc/))2 • 

More in general we can see that the holonomy operators fll27p with k = 0, that is 
i^Pi,ki,mi;0;p2,k2,m2 generate a closed subspace. Hence the dynamics only involves this 
subspace of the universal boundary Hilbert space. Indeed the integrand of the partition 
function, that is the effective face weights do show an enhanced symmetry in the case 
of the Barrett Crane model: The ^'-function is given as E{g^,g^) = ^ su {2){9^ {9^)~^) 
where {g^,g^) G SU{2) x SU{2). It is invariant under SU{2) multiplication (of the 
diagonal subgroup in SU{2) x SU (2)) from the left and from the right. This is a stronger 
symmetry than in the general case where E is just required to be invariant under the 
adjoint action of the subgroup H. Because of this enhanced symmetry the effective face 
weights just depend on variables G x G/H, associated to every half edge. The subspace 
of the Hilbert space spanned by basis states with /c = is exactly the subspace invariant 
under this H group action. 

Let us turn to the W operator and see in which sense it is a constraint on the curvature 
(as is the case in BF theory). The effective face weights can be computed to 



Uj'f 




N 



^H{Y[ljHejVjHvjej+,l/ ) ■ (156) 



J=l 

Here 7/ are group variables in H = SU{2) and abbreviates gye{9ve)~^ ^^e product of 
left and (inverse) right copy of the group. N denotes the number of edges in the face /. 

The requirement for a non-zero face weight is that there exists a set of group elements 
k'j, which have to be in the conjugacy class of := Heivri^viej+i, such that the product 
of the k'j is equal to the identity. For faces with more than two edges this condition is 
generically satisfied and the effective face weights only vanish on measure zero subman- 
ifolds in the configuration space. Thus we cannot read off a curvature constraint for the 
general case. 



6.6 Discussion of transfer operator and its limit 

A question of intense research is the relation between the dynamics as defined by the 
spin foam models on the one hand and the dynamics as defined by the Hamiltonian 
constraints in loop quantum gravity [98] on the other hand [771 [H El EH [37] . The 
model example is 3D gravity, which is equivalent to BF theory [77]. In this case the 
transfer operator is a projector and its image can be described by quantum constraints 
which reflect the diffeomorphism symmetry of the model. These constraints can also be 
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encoded into recursion relations (in the spin representation) which can be interpreted as 
the Wheeler DeWitt equations of 3D quantum gravity [22l [38] . 

In 4D gravitational models the situation is complicated by two main issues. One 
is that diffeomorphism symmetry in discrete 4D gravity models are generically broken 
even on the classical level, so that one cannot expect the transfer operators to be a pure 
projector and to lead to constraints [lH [131 SB]- There is an exception to this general 
picture, in cases that the dynamics allows only for fiat geometries. This is the case for 
special triangulations described in \5T\ HH] and more specifically for so called tent moves 
at four-valent vertices [IS], on which we will comment more below. 

The second main issue is the appearance of simplicity constraints. The natural bound- 
ary Hilbert spaces for the current spin foam models are based on G = 5*0(4) (or 
G = 5*0(3,1)) holonomies, whereas the LQG Hilbert space is based on H = SU{2) 
holonomies. Although we can formally define the transfer operator (i.e. for the EPRL 
model) as an operator on the LQG Hilbert space it rather involves the multiplication 
operator W with G holonomies. 

Let us comment more on the two issues and point out avenues for further research. 
In case of broken diffeomorphism symmetry the transfer operator is rather an evolution 
operator corresponding to some finite time step instead of a projection operator. Thus 
one could ask for the limit in which this finite time step is taken to be small in order to 
extract the time generator, that is the Hamiltonian. Indeed this is the standard proce- 
dure for lattice gauge theory, in which the time step is encoded in the lattice constant 
(in time direction). In gravity however we do not have such an explicit lattice constant 
(or other coupling constant). Rather the time distance is boundary state dependent and 
might only emerge semiclassically. That is to take the limit of infinitesimal time steps 
we need to entangle this procedure with some semi-classical limit. We discussed different 
boundary Hilbert spaces, hence different types of semi-classical states are possible, in- 
cluding Hall-like semiclassical states [99l [19] adapted to the universal boundary Hilbert 
space and coherent states used in the discussion of the semi-classical limit of spin foam 
models [2H1 12S1 [23 [271 [Ml m]. 

An alternative to considering the limit of infinitesimal small time like distances is to 
'improve' the transfer operator by basically coarse graining [lll[I71[l3], i.e. considering 
an effective transfer operator T' = . This again is similar to standard statistical 
lattice theories, where the limit for — oo leads to a projector on the eigenspace 
corresponding to the highest eigenvalue of T. However, in general some scaling is re- 
quired in order to obtain an interesting projector, i.e. the highest eigenvalue should be 
rather highly degenerate. To this end it might be necessary to refine also the spatial 
discreteness, as the limit of continuous time but discrete space might rather not lead to 
a restoration of diffeomorphism symmetry |75] . 

Here we discussed a global transfer operator which acts on the entire hypersurface. 
An alternative, more adapted to the 'multi-fingered' time evolutions are so-called tent 
moves, discussed in [961 HH [131 [ISj, which rather evolve single vertices in the hypersur- 
face. To define the corresponding transfer operator one needs to provide the possibility 
of gluing wedges, which is discussed in section [H Tent moves are especially interesting 
as these would lead to localized, vertex based (Hamiltonian) constraints (in case the 
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symmetries are realized). Indeed for classical Regge calculus tent moves at four-valent 
edges lead to constraints as the simplicial geometery remains flat in this case [HI |19] 
and diffeomorphism symmetry is preserved. Here it would be interesting to know if the 
spin foam transfer operator corresponding to a tent move at a four-valent vertex leads 
to a projector and therefore constraints or not. More generally the question is whether 
we can obtain a dynamics describing fiat geometries for triangulations that in Regge 
calculus only allow for fiat metrics due to topological reasons [5T] . 

The second more technical main point concerns the simplicity constraints. As we have 
seen the transfer operator naturally involves as a G-group holonomy and boundary 
Hilbert spaces based on G. Hence, even if we can formally define matrix elements of the 
transfer operator on the LQG Hilbert space (based on H holonomies), a comparison to 
a (Hamiltonian) operator expressed in terms of H holonomies and fluxes is rather diffi- 
cult. An alternative is provided by the recent work [331 El] which provides a canonical 
connection formulation based on G holonomies. In this case the Hamiltonian constraints 
have to be augmented by a term that makes them gauge invariant with respect to the 
(primary) simplicity constraints. One could argue that this term is taken care off in spin 

1 /2 

foams, as (or Ke or /i) projects back onto the solutions of the (primary) simplicity 
constraints. However, heuristically the transfer operator is the exponential of the Hamil- 
tonians. The additional term in the continuum Hamiltonian constraints takes care of 
staying on the simplicity constraint hypersurface at all times. In contrast, with the spin 
foam transfer operator we rather project onto the simplicity constraint hypersurface in- 
between discrete time steps. The discrete time steps itself involve G group holonomies 
which in general map out of the subspace defined by the simplicity constraints. This 
is reminiscent of discussions in [1], which argues that simplicity projectors should be 
inserted at each point of the G-holonomies. More generally this problem is connected 
with the issue of how secondary simplicity constraints [S], IIUUI \5T\ \5U\ are implemented 
into spin foams. The work here offers the possibility to check in which sense the impo- 
sition of primary simplicity constraints at consecutive time steps leads to an imposition 
of secondary simplicity constraints, as for instance argued in [73]. These insights might 
allow a relation between G holonomies and holonomies involving the Ashtekar-Barbero 
connection [HZ] , which in turn will ease the comparison of the transfer operator with the 
Hamiltonian constraints. 

7 Discussion 

We have shown that the new holonomy representation presented in this and the compan- 
ion paper [15] provides several advantages. There is a clear parametrization of the space 
of models. From the parameters one can easily read off the reality of the amplitudes, 
and it turns out that the BC, EPRL, FK and BO model have real amplitudes in this 
representation. The holonomy models lead to a natural boundary Hilbert space, which is 
introduced in this work. This universal boundary Hilbert space is the same Hilbert space 
for different choices of simplicity functions E. Thus different possibilities to impose the 
simplicity constraints can be compared in one and the same Hilbert space. Furthermore, 
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as the models can be naturally defined on arbitrary 2-complexes we also can obtain ar- 
bitrary (and not only 4-valent) graphs on which the boundary Hilbert spaces are based. 
A (dynamical) notion of cylindrical consistency, related to coarse graining, can also be 
introduced on these Hilbert spaces [121 SS] . 

The universal boundary Hilbert space results from a definition of 'equal time' slices 
which is most natural for lattice gauge theory. Alternatively we can adopt slices and 
gluings more custom to spinfoams and obtain as boundary Hilbert space the Hilbert 
space of projected spin networks. We detailed the conditions under which such map be- 
tween the different slicings is possible and constructed the corresponding /i-map. These 
considerations allowed also a discussion on the different kind of basic building blocks, 
i.e. faces, half faces and wedges, and the possible gluings which allow a combination of 
these building blocks to the full partition function. 

We explicitly constructed the representation of different current spin foam models 
in our new representation - the main difference between the models are the simplicity 
functions E, that imposes the simplicity constraints. 

Finally, we derived a general form of the transfer operator on the different boundary 
Hilbert spaces. On the universal Hilbert space this form is given by T = K^^^ PqW PgKq^ ^ 
with W representing a product over holonomy operators over faces, Pq is the projector 

1 /2 

on the G invariant subspace and Kq can be seen as imposing the simplicity constraints 
on the Hilbert space. Similarly on the Hilbert space of projected spin network func- 
tions we have T' = fi^Wfi with taking over the role of K to impose the simplicity 
constraints. This form might shed some light on the discussion of how to best impose 
simplicity constraints into spin foams O |50l |3ll |61] . In the current models the secondary 
simplicity constraints are not imposed. It is argued that the imposition of primary sim- 
plicity constraints on each time step should also lead to the imposition of secondary 
constraints. (Note however that even the primary constraints are imposed weakly in 
the EPRL type models.) Indeed, here this is made obvious in the form of the transfer 
operator, who is projected by either Kq or /i. Of course the projections are only in- 
serted at discrete time steps and not continuously in time. Related to this issue is that 
fundamentally the transfer operator still includes a, G = S'0(4)-holonomy operator and 
not a holonomy operator based on the Ashtekar-Barbero connection. A question for 
future research is how this form relates to a discretion which starts from the beginning 
with the 5'f/(2)-Ashtekar Barbero connection, which arises by classicaly solving for the 
simplicity constraints [20] • This question can be also studied in a semi-classical limit 
using semi-classical states in either of the boundary Hilbert spaces or the techniques 
presented in [65] . 
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